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In this paper, we discuss the role of an external magnetic field on the dy- 
namically generated fermion mass in even-flavor QED in three space-time 
dimensions. Based on some reasonable approximations, we present analytic 
arguments on the fact that, for weak fields, the magnetically-induced mass 
increases quadratically with increasing field, while at strong fields one crosses 
over to a mass scaling logarithmically with the external field. We also confirm 
this type of scaling behavior through quenched lattice calculations using the 
non-compact version for the gauge field. Both the zero and finite tempera- 
ture cases are examined. A preliminary study of the fermion condensate in 
the presence of magnetic flux tubes on the lattice is also included. 



I. INTRODUCTION 

The particle mass generation via dynamical symmetry breaking has been a much-studied 
scenario in particle physics as well as in condensed-matter systems. In the recent years this 
phenomenon has been studied in the presence of background fields, such as constant external 



magnetic fields p]]- ||, following and extending the formalism developed by Schwinger [ITO 
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The formalism has been applied to models that had gauge and/or four-fermion interactions. 
It was found that such constant background configurations can enhance the dynamical sym- 
metry breaking by driving the critical coupling to a smaller value and thus catalyzing the 
symmetry breaking. A concrete example of this phenomenon, of relevance to us in this 
work, is the dynamical chiral symmetry breaking of chiral symmetry in massless QED ( in 
three and four dimensions) in the presence of an external magnetic field where the 

dynamically generated fermion mass depends on the value of the external field. 

The magnetically catalyzed mass generation in (2+1) dimensional QED may have in- 
teresting condensed-matter applications U\M, given the suggestions that high-temperature 



superconductors can be described effectively by field theories like QED 3 [|TTJ or by non- 
Abelian gauge models based on the group SU(2) x U(l) fT2l,|7H []• Indeed, there is experi- 
mental evidence for the opening of a second (superconducting) gap at the nodes of the gap 



in certain d-wave superconductors in the presence of strong external magnetic fields |T3 
As remarked in |TJ], in the context of condensed- matter- inspired models, the scaling of the 
thermal conductivity with the external field is different between the gauge [|7],[| and four- 
fermion theories (T^|. Thus, a detailed study of the magnetically-induced chiral symmetry 
breaking phenomenon in the context of QED3 is phenomenologically desirable, given that 
such studies may lead to more detailed experiments in the spirit of JIB] , that can probe deep 
in the structure of the novel high-temperature superconductors. 

In (2+1) dimensions chiral symmetry can be defined only if the number of fermion flavors 
is even This fact is relevant for a planar high-T c superconducting antiferromagnetic 

system flTT| , |l7|l which comprises of two sublattices. Within a generalized spin-charge 
separation framework |Tj|, there will be two species of charged fermion excitations (called 
holons), one associated with each sublattice [11, Finally, the (2+1) dimensional theory 
with even number of fermion flavors |7| can be viewed as a dimensional reduction of the 
four-dimensional effective Lagrangian of ||. 

In QED3, the magnetic catalysis of the chiral symmetry breaking for strong external 
fields is established by looking at the Schwinger- Dyson equations ||[7|]- In these works 
the Landau level formalism was used to truncate the fermion propagators to the lowest 
Landau level. This formalism is satisfactory for certain aspects of the magnetic catalysis for 
strong magnetic fields J7j, but for weak fields the result can definitely be questioned, given 
that in that case the spacing between Landau levels becomes small, and one effectively 
deviates from the lowest Landau level description. Recently two of us [|19[] have looked at 
the role of higher Landau levels and showed that they contribute by inducing a (parity- 
violating) magnetic moment which scales with the applied magnetic field. Moreover the 
role of higher Landau levels in inducing a critical temperature even in the free fermion 
case, under certain circumstances, was emphasized in ||. For all of the above reasons it is 
important to incorporate the effects of all the higher Landau levels in the Schwinger- Dyson 
formalism, avoiding the use of the mean field Landau level decomposition altogether. This 
is what we shall attempt to do in the first (analytic) part of this paper. We shall compare 



1 The relativistic (Dirac) nature of the fermion fields is justified by the fact that they describe the 
excitations about the nodes of a d-wave superconducting gap. 
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these results by performing some preliminary (quenched) lattice analyses in the second part 
of the paper. In the latter respect we have to mention that the quenched approximation for 
fermions employed here allows for the ladder gauge quantum fluctuations in the fermion free 
energy to be incorporated, but prevents the use of internal fermion loops, as the treatment 
of the latter requires an algorithm for treating dynamical fermions which is currently under 
construction. 

The paper is organized as follows. In section 2 we give a brief review of the SU(2) x Us(l) 



model of |T2| , as well as the Dirac algebra in three-dimensional spacetime with an even 
number of fermion flavors. In section 3 we review the Schwinger-Dyson (SD) equation for 
the fermion propagator in the absence of the the external magnetic field. In section 4 we 
present the results for the case of strong external fields, where a logarithmic scaling of the 
induced condensate with the external magnetic field occurs. In section 5 we present the 
SD equations for the weak magnetic fields ignoring the photon polarization to make contact 
with the lattice result presented in the second half of the paper. We show that under certain 
approximations, the scaling behaviour of the condensate with the external magnetic field 
can be found. In the next section, we attempt to go beyond the quenched approximation 
analytically, by including the photon polarization and modify accordingly the Schwinger- 
Dyson equations. The analysis becomes very complicated to be handled analytically for 
finite temperatures, and this is the reason why we turn to the lattice formulation of the 
problem in section 7, where we set up the formalism and relevant notations. In section 8 the 
lattice results are presented for both zero and finite temperatures; in addition, a preliminary 
extension of the results to the non-uniform magnetic field cases is attempted by examining 
the magnetic catalysis phenomenon in the case of flux tubes. This model may constitute 
a prototype for the study of the effects of electromagnetic vortices in condensed matter 
systems, which are of relevance to high-temperature superconductivity. Conclusions and 
outlook are presented in section 9. 



II. THE MODEL AND ITS SYMMETRIES 



The SU(2) x U(l) model of |2(J is a toy model for dynamical electroweak gauge symmetry 
breaking in three dimensions, while in the context of condensed-matter systems, the SU(2) x 



Us (I) model of [12| is based on a gauged particle-hole symmetry, via a suitable extension 
of the spin-charge separation The holons transform as a doublet under the SU(2) 

(particle-hole) symmetry. In this respect the model is different from other SU(2) x £7(1) 
spin-charge separated theories, which are based on either direct gauging of genuine spin 
rotation SU(2) symmetries PTf , or non-Abelian bosonization techniques P^J2^| . The phase 
diagram of the model of [Q, and the associated symmetry-breaking patterns, are quite 
different from these other models. 

The three-dimensional continuum Lagrangian of the model is given ( in Euclidean metric, 
which we use hereafter) by [F2U|.|7[. 



£ = -i(i^) 2 - -{Q^Y + VtD^i - m^^i (2.1) 

where D M = <9 M — igia^ — ig 2 a a B a ^, and F^, Q^ v are the corresponding field strengths for 
an abelian ("statistical") Us(l) gauge field a s and a non-abelian ("spin") SU(2) gauge field 
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B®, respectively. Due to the antiferromagnetic nature of the condensed matter system the 
fermions are four-component spinors, i = 1, • • • N. We note that may be written as 



*< - V' ■ (2-2) 




Then the Lagrangian decomposes into two parts, one for and one for \I/j2, which will 
be called "fermion species" in the sequel. The presence of the even number of fermion 



species allows us to define chiral symmetry and parity in three dimensions [IB|, which we 
discuss below. The bare mass m term is parity conserving and has been added by hand 
in the Lagrangian ( |2.1| ). In the model of ]T2]|7|, this term is generated dynamically via the 
formation of the fermion condensate < > by the strong Us{l) coupling. However, for 
our purposes, the details of the dynamical mass generation is not important and hence it 
will be sufficient to include a bare mass term for the holons representing the mass generated 
by the (strongly coupled) Us (I) interactions in the superconducting phase. 

In what follows we shall ignore for simplicity the non abelian gauge group structure 
and concentrate only in the Abelian model in the presence of an external electromagnetic 
field, which should not be confused with the statistical abelian gauge field Us(l)- The 



incorporation of the gauged SU(2) structure leads to a much richer phase structure [|24.19 
and we reserve the discussion for future publication. 

For even-flavour models a convenient representation for the 7 M , fi = 0, 1, 2, matrices is 
the reducible 4x4 representation of the Dirac algebra in three dimensions |TB 





(2.3) 

where a are 2x2 Pauli matrices and the (continuum) space-time is taken to have Euclidean 
signature. 



As well known [16| there exist two 4x4 matrices which anticommute with 7 M , /i = 0, 1, 2: 
73= (° I), * = \) (2.4) 



where the substructures are 2x2 matrices. These are the generators of the 'chiral' symmetry 
for the massless-fermion theory: 

* -> ezp (26*73)* 

* -> exp(iuj 5 )^. (2.5) 

Note that these transformations do not exist in the fundamental two-component represen- 
tation of the three-dimensional Dirac algebra, and therefore the above symmetry is valid for 
theories with even fermion species only. 

For later use we note that the Dirac algebra in (2 + 1) dimensions satisfy the identity: 
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yy = - r 3 e^ A 7 A ; r 3 = z 7 37s = ( q 1 _^ j 

'T ^ 'T — — 'T 

rfrfrfrf = _ 7-3/^' 

yy-yyy = _^ 7 fc _ ^y + <p*y (2.6) 

which is specific to three dimensions only. Here the Greek indices are space time indices, 
and repeated indices denote summation. 

Parity in this formalism is defined as the transformation: 

P: ^(xW 2 ) -> ^^(i ,-! 1 ,! 2 ) (2.7) 

and it is easy to see that a parity-invariant mass term for \l/ amounts to masses with opposite 
signs between the two species |l(J, while a parity- violating one corresponds to masses of equal 
signs. 

The set of generators 

Q = {1, 73, 75, A = z 7375 } (2.8) 

form ^U|,|T2J a global £7(2) ~ SU(2) x £7s(l) symmetry. The identity matrix 1 generates 
the Us(l) subgroup, while the other three form the SU(2) part of the group. The currents 
corresponding to the above transformations are: 

Jl = ^ 7At r^ T = 7 3, 75, ^7375 (2.9) 

and are conserved in the absence of a fermionic mass term. It can be readily verified that 
the corresponding charges Qr = J a^x^/^T^ lead to an 577(2) algebra p0| |: 



[Q 3 ,Q 5 ]=2iQ A [Q 5 ,Q A ] = 2iQ 3 

[Qa,Q 3 ] = 2zQ 5 (2.10) 

In the presence of a mass term, these currents are not conserved: 

= 2m*r* , (2.11) 

while the current corresponding to the generator 1 is always conserved, even in the pres- 
ence of a fermion mass. The situation is parallel to the treatment of the SU(2) x SU(2) 
chiral symmetry breaking in low-energy QCD and the partial conservation of axial current 
( PCAC). The bilinears 

Ai = ^73^, A 2 = ^75^, ^3 = 

= *7m73*, B 2ll = ^75*, B 3ft = W 7m A#, // = 0, 1, 2 (2.12) 

transform as triplets under SU(2). The SU(2) singlets are 

A = *A^r, B i>fl = Vj^ (2.13) 
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i.e. the singlets are the parity violating mass term, and the four-component fermion number. 

We now notice that in the case where the fermion condensate A3 is generated dynamically, 
energetics prohibits the generation of a parity-violating gauge invariant SU(2) term P5] 



and so a parity-conserving mass term necessarily breaks [|7| the SU(2) group down to a 
T3-C/(l) sector ]TI[, generated by the 03 Pauli matrix in two-component notation. Upon 
coupling the system to external electromagnetic potentials, this phase with massive fermions 
shows superconductivity. The superconductivity is strongly type II pT|J7[| as the Meissner 
penetration depth of external magnetic fields turn out to be very large,^ and hence the study 
of the response of the system to the external electromagnetic fields is justified. 

III. THE SCHWINGER-DYSON EQUATION FOR THE FERMION 

QED 3 is a superrenormalizable theory which is confining in the infrared regime. Ac- 
cordingly, it acts as a simple prototype for the analysis of the chiral symmetry breaking in 
QCD. The standard tool for investigating the chiral symmetry breaking are the celebrated 
Schwinger-Dyson equations. In this section, let us set up the Schwinger-Dyson equations 
for the fermion propagator. 

The Schwinger-Dyson equation concerning the fermion propagator Sf{p) (for zero bare 
fermion mass) is given by: 

d 3 k 

where T u is the fermion-photon vertex function and D^ v is the exact photon propagator. 
However, to this order, let us make the following approximations: 



f d k 

s- F \v)=i-v-g j j^— 3 rS F (k)T^k,p-k)D, u (p-k) (3.1) 



1. Use the bare vertex function, namely 

T v (k,p-k)= gi \ (3.2) 

so that the gap equation reads: 

/d J? 
^S F (k)Y(k,p- k)D, v {p - k) (3.3) 

2. Now, we choose the following ansatz for the full fermion propagator: 

S~\p) = A(p) 7 % + B(p)j ■ p + E(p) (3.4) 

Using this ansatz, let us now perform a trace over the gamma matrices in ( |3.1| ). This 
gives us the following gap equation 

= 9 J J^f A-kl + + S'(fc) ^ D ^ P ~ k) (3 - 5) 



The high-temperature superconducting oxides are strongly type II superconductors. 
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3. To further simplify the gap equation let us use the zeroth order result for the wave- 
function renormalization, namely A(p) = B(p) = 1, which is often justified in the large 
iV argument |16|] (see however [PT| ) so that eq.(|3.5|) reads: 

4. The photon propagator D fiu (k) can be replaced by the ladder resummed propagator 
which can be justified in the large- N limit. The resummed propagator (in the absence 
of the magnetic field) is given by 

n ( n ) - eL_ - \JZ e!_ (q 7^ 

^w-p2 (1 _n(p))- p2 (1 + g) C3 " 7; 

The gap equation thus obtained in the absence of the magnetic field can be solved using 



the bifurcation method [16]. There are two solutions namely, 

^(p)-p-^ 2 \ MP)~P 1 -*^ N (3-8) 

where N is the (large) number of fermion flavours. However, it is natural to expect that 
these solutions will change in the presence of the external magnetic field; we will discuss this 
generalization in the following sections. 



IV. THE DYNAMICALLY GENERATED FERMION MASS AT STRONG 

MAGNETIC FIELDS 

As mentioned above, hereafter we consider only the abelian gauge group Us(l) in the 
presence of an external electromagnetic potential A e ^ f , corresponding to a constant magnetic 
field B, perpendicular to the spatial plane. The dynamics is described by the Lagrangian: 

C = -\{F^f + tfZWP - m*tf (4.1) 

where = d^—iga^—ieA^. The mass m here should be viewed as an (infrared ) regulator 
mass. In the dynamical mass generation scenario investigated below via the SD method m 
should be set to zero, given that the dynamics of the gauge field and the magnetic field are 
both responsible for the appearance of a mass in the fermion propagator. For the lattice 
analysis, on the other hand, the presence of an initial small 'bare' regulating mass m ^ 
appears necessary ||. 

We commence our analysis by noting that the presence of an external magnetic field, 
perpendicular to the spatial plane X1X2, breaks Lorentz and translational invar iance. The 
configuration space form of the fermion two-point function G(x, y) for the three-dimensional 
problem at hand has the generic form ||: 

G(x, y) = exp (|(x - y)Mf (* + y)) G(x - y) (4.2) 
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where A e ^ denotes the external electromagnetic potential, corresponding to a constant ho- 
mogeneous magnetic field perpendicular to the spatial plane Xix 2 : A^ = (o, — -|x 2 , ^X\ 
(in an obvious notation). The field-dependent phase factor in ( f4.2|) breaks translational in- 
variance, implying that, in general, G(x,y) does not admit a Fourier transform expressible 
in terms of a single momentum (vector) variable k. 

The translationally-invariant part G(x—y) has a Fourier transform Sp(k) of the form [ITU 



S F (k) = i [°° dse- s(k o +m2+k 
Jo 



2 tanh z 



) [(m — 7 • k) — z(7i&2 — 72&1) tanh z] (1 — 27172 tanh z) 

(4.3) 



where m is the mass of the fermion, and z = s eB. Note that we are distinguishing between 
the coupling constant g for the statistical U(l) gauge field and the electromagnetic charge e. 



The Schwinger propagator admits the following expansion in terms of the Landau levels |26 

:-i) n D n (k ,k) 



S F (k) 



te eB 



E 



q kl + m 2 + 2enB ' 



(4.4) 



with 



D n (k Q ,k) = 
(m - 7qA;o) 

(m - 70/co) 



,2k 2 
~e~B 



n-l\ 



,2k 2 
~eB 



2k 2 
~~eB 



+ ir 3 7o L n ( 



~eB 



,2tf_ 
~eB 



+ 4( 7 • k)L\_ x 
+ 4( 7 .fe)4_ 1 ( 



2fc 2 
^~e~B 



(4.5) 
(4.6) 



For QED 3 the scaling of the dynamically generated fermion mass with the external 
magnetic field had been discussed by Shpagin || and two of us J7J. Let us begin with 
the case when the external magnetic field is very strong. As stated in the introduction, in 
this case it is sufficient to truncate the fermion propagator ( in the absence of the Us (I) 
interactions) to the lowest Landau level ( |4.4j ), so that we get: 



qLLL 
Dp 



-hi 

ie eS 



1 



m + 7 k° 







27 1 7 2 



(4.7) 



As we will be dealing with the lowest Landau levels only, it is expedient to choose the ansatz 
for the "exact" propagator for the lowest Landau level fermions to be of the form: 



qLLL 
Dp 



_ hi 

ie <= B • 



1 



Hence, following 



E(fc) + A(k)^k® 
the gap equation for the lowest Landau level fermion is given by 



d 3 k 

W) 



E(fc) 



A 2 kl + Y?(k) 



D 00 (p-k) 



(4.8) 



(4.9) 



According to |J the photon vacuum polarization gets suppressed as ^7= at strong mag- 
netic fields and the photons become almost free. Thus in the strong field limit the photon 
propagator is given in the Landau gauge by the expression: 



8 



D Av) = ^ p ^ ( , \ ( 4 - 10 ) 



P 2 p 2 (l + 0.14037- VejI> 
Accordingly, we have 



w x ~ 2 f j2 U -4 1 fi (Po- k o) 2 \ f dk o (A i i \ 

E (p) = » £B i 1 ~ J J A*k% + ' (411) 



where g 2 = — — Let us set p to zero. Then, 



1+0.14037-2=— 



„. . f k 2 , 9 d 2 k r (i/Cn £(&) 

s (°)=»/ e " 3k MTw/(4^riw (412) 

For strong fields e5 — > oo, we suppose that setting ~ S(0) and A ~ 1 yield a sufficiently 
good approximation p7 |. Setting k 2 = x, the gap equation becomes: 

1 = hS dk °S {kl I xf kl + E 2 (0) " (413) 

Assuming that the dynamically generated fermion mass is much smaller than the external 
magnetic field, i.e. S(0) << \J eB, we cut off the x integration by y/ eB and after the k 
integration we get the transcendental equation: 



m „£j^te-i(i-?m + m) ,4,4) 

Air je(o) \y y z V 



The final result reads: 



S(0) - 2a In (^f ) +0(^§), where a = £. (4.15) 



This equation can be solved numerically as in J/J, leading to a logarithmic scaling of the 
induced fermion condensate with the magnetic field, S(0) ~ ln(V eS / a) . 

Note that the most important aspect of this type of behavior comes from the presence 
of the exponential factor in the form of the free propagator in (|4.7|) . However, when the 
external magnetic field is weak, one has to include all the higher Landau levels as the levels 
become closely spaced. Also, the wavefunctions for these levels grow with momentum as they 
involve Laguerre Polynomials. Hence, one has to work with a generic ansatz for the fermion 
propagator, such as ( |3.4j ). However, in view of the breaking of translational invariance by 
the field- dependent phase factor in (|4.2|) , a straightforward application of this ansatz is not 
possible. Nevertheless, as we shall discuss below, such an ansatz can still give qualitatively 
correct predictions for the scaling of the induced condensate with the external field. 
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V. THE DYNAMICAL FERMION MASS IN WEAK EXTERNAL MAGNETIC 
FIELDS UNDER QUENCHED APPROXIMATION 

We are looking at the leading scaling behaviour with the magnetic field intensity of the 
dressed fermion propagator in the presence of an external magnetic field eB for the case 
of weak fields eB << £(0) where £(0) = m is a dynamically generated mass due to the 
statistical Us(l) interactions in the model. Obviously, since m oc g 2 , whereas g denotes the 
coupling of these interactions, the weak field limit is achieved for relatively strong gauge 
interactions. However, as stated earlier, we are interested in the behavior of the system 
under weak magnetic fields as well and in this regime the Landau level decomposition is not 
particularly helpful. 

This is a technically involved problem, and we shall not attempt to solve it exactly in 
what follows. Instead we shall make an attempt to present an approximate treatment, which 
hopefully captures the important qualitative features of the phenomenon. As we shall see 
later on, preliminary lattice calculations will support the results obtained in this section. 

The main complication arises from the fact, mentioned in the previous section, that 
the presence of an external field, perpendicular to the spatial plane, breaks, in addition 
to Lorentz invariance, also translational invariance in the spatial plane. This is apparent 
from the configuration space form ( |4.2| ) of the fermion Green's function. The breaking of 
translational invariance is manifested through the field-dependent phase factor in ( |4.2|) . A 
Schwinger-Dyson equation for G(x,y) can then be obtained as in ||, but unfortunately, 
due to the form ( |4.2| ) , a passage to momentum space with an appropriate Fourier transform 
based on a single momentum variable is not feasible. 

Below we shall make a modest attempt to calculate analytically the scaling behaviour 
of the chiral condensate with the magnetic field, at least qualitatively, and then compare 
with the lattice results. This is possible by adopting an ansatz for the fermion propagator 
in momentum space, S^fjo), as if the translational invariance breaking phase factor in ( |4.2|) 
was absent. Specifically, we assume that the fermion propagator in the presence of a weak 
external magnetic field can still have a Fourier transform, based on a single momentum 
variable, which takes the approximate form: 



All the qualitative information about the effects of the external field B is encoded inside the 
coefficients. As we shall see later on, by comparison with the preliminary lattice results in the 
quenched approximation, the qualitative features of the scaling of the magnetically-induced 
chiral condensate with the external field seem to be correctly captured by the above ansatz. 
It is understood, of course, that these results should be taken with caution and viewed only 
as preliminary. A complete lattice analysis involving dynamical fermions, which will settle 
these issues, is still pending. Hopefully it will constitute the topic of a forthcoming work. 

Under the above approximation, the Schwinger-Dyson equation for the fermion propa- 
gator in the massless limit assumes the form: 



s f 1 (p) = l A o + -4i7i72]7oPo + C 7 • p + [S + £17172]. 



(5.1) 




(5.2) 
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It is sufficient for our purposes to use the bare photon vertex (T u = j u ), and set the 
wave-function renormalization to one. These will be justified later on. 
Taking the trace in the above equation, we obtain: 

4S(p) = -g 2 J 0r 3 D,u(p - k)Tr (^S F (k)Y) (5.3) 

where we denoted E by S(p). A more important remark is that inside the integral we have 
approximated the (unknown) fermion propagator S F (k) by the form S F (k) of equation ( [4. 3D , 
which is the propagator in a homogeneous external magnetic field. One has: 

Tr (YS F {k)Y) = J™ rfse- s(fc « +fc2 ^ +E2) Tr[ 7 ^(-7 • k + T.{k))-f + 

ry M (— + 7 ■ /c)7i727 I/ tanh,2 — ry M ( 7 ifc2 — 7 2&i)tanh;z(i — i 7 i 72 tanh2;) 7 iy ] (5.4) 



Since we work in the Landau gauge, D^q) is given by q d ^ 4 "^\ Using the following 
identities of the 7 matrices in Euclidean space: 

Tr(jrf) = -12 

(q ■ if = -<f 

Tr^O^D^q) = ~\^0 (5.5) 
where O is any operator, one can write ( |5.4| ) in the form: 

D, v (p - k)Tr (YMk)Y) = -T^p f e- s(fc ° +fc2 ^ +E2(fc)) . (5.6) 
Using (|5.6j), ( p.3[ ) becomes: 

^ = W?l« dS l J^me-^'^^ (5.7, 
Setting p = and approximating ~ S(0) = m in the integrand yields: 

292 r ds [d 3 k e -.(*3+*»*^»-H^) > ( 5. 8) 



(27r) 3 7o 7 A^ + fc 2 
Now we use the parametrization: fco = kcosd, k 2 = k 2 sin 2 6, and get: 

fp- poo r+1 g—sm? 

1 = — 3- / ds dx . , (5.9) 

4:712 Jo J-i s ( x 2 + (i_ x 2\t™hz\ 



with x = cosO. Then the x integration can be performed: 

1 g 2 f°° , _™2 1 v/i + ^Jz - tanhz 

' cfee z e B — log-^-p , (5.10) 



47r5 v r eB Jo \J z — tanhz yfz — yfz — tanhz ' 
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where we have set z = eBs. At this point we define two dimensionless variables: // 



2tt%, / = 2n^r, and write the last equation in terms of /i and / 
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1 = ~ fdze-^ , 1 fog ^+V^ggf , (5 . n) 
2^n f Jo \/ z — tanhz y'z — y z — tanhz 

For each value of / the above relation gives a corresponding value for /j. The result is a 
universal curve //(/), which we depict in the figure |3] as a dashed line. This result can be 
translated very easily to the language of dimensionful parameters, as well as to the lattice 
parameters. Notice that up to now there is no restriction to weak fields. 

One may also derive an analytical approximation for the regime of weak magnetic fields. 
Starting from eq. Q5.8Q we write the term e~ sk in the form e~ sk ■ e~ sk (- SI r £ ~ 1 \ expand 
the second exponential in a power series of z = eBs for weak fields and retain the terms 
which are at most of sixth order in eB. Then the same parametrization as before is used 
and the integrations over k and x are similarly carried out. Then the equation involving the 
integral over z is replaced by: 



g 2 

l = — 

2tt 



1 1 e 2 B 2 7 e i B A 773 e 6 B 6 ^ . s - 

_i_ _i_ _i_ 

m 12 m 5 48 m 9 960 m 13 1 



(5.12) 



For weak fields we may solve equation fl5.12j ) by substituting the expansion m = m + 
rri2 (eB) 2 + m 4 (e-B) 4 + ... in (|5.12p and equating the coefficients of equal powers of eB to 
determine the coefficients m&. The resulting solution of (|5.12 ) is: 



g 2 ( 4n A e 2 B 2 400n 8 e A B A 534208vr 12 e 6 J B 6 ^ /8r , H A 

m=— 1+ 5 r R + 53 + 0(e 8 B 8 ) , (5.13) 

2tt V 3^ 8 9g w 135g 24 J 

which may also be written as: 

"- 1 + ^ 4 -^ 8 + ^ 12+ w 16 )- <» 4 » 

The above relations show that for weak magnetic fields the dynamically generated mass is 
quadratic in B. For somewhat bigger magnetic fields, however, the quadratic behaviour is 
compensated by a negative quartic contribution and the increase with the magnetic field 
resembles very closely a linear dependence. Of course for even bigger magnetic fields higher 
order contributions take over. 

Although we have already depicted in figure ^ the full solution //(/), we also plot the 
solution of equation ( |5.12j ) (we note here that we did not actually use the solution ( |5.13|) , 



but rather the one which results if (|5.12|) is quenched to order B 2 before solving); this is 



done to gain some feeling about the accuracy of the quadratic approximation, which will be 
the only possible approach in the case where also dynamical fermions are taken into account. 
The approximate solution is the full line and is restricted in the region of small B, where it 
is appropriate; it is quite good up to y^pr ~ 0.1. 

As we have just seen, in the quenched approximation one uses the free photon propagator 
as the fermion loops ( which modifies the fermion propagator ) are ignored. In the next 
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section, we examine whether the inclusion of the photon polarization modifies the scaling 
behaviour of the gap function with eB discussed above. 

Before closing this section two important remarks are in order. First, it should be noted 
that the presence of the translational invariance breaking phase factors in Q4.2|) , which have 
been ignored in the above treatment, will affect the numerical coefficients of the even powers 
of B in (|5.13|) . This can be seen easily from the form (}4.2|) by an expansion in powers of 
the weak field B (and restriction to the real part, assuming hermiticity of the translational 
invariant parts). The important issue is the sign of the various terms. As we shall see 
later on, comparison with the (quenched) lattice results confirms the scaling with B given 
in ( |5.13| ), thereby justifying the above approximate method of dealing with the problem, 
at least for qualitative purposes. Second, the above analytic treatment, leading to ( |5.13D , 
was based on the approximation of replacing S(/c) inside the integrals in the pertinent 
integral form of the gap equation by a constant S(0) = m, the so-called "constant-mass 
approximation" . This is sufficient for the qualitative purposes of this work, where the main 
interest lies on the scaling of the induced condensate with the magnetic field. It should be 
remarked though, that attempts to go beyond the "constant mass approximation" have been 
made in the literature, specifically in the context of three-dimensional QED in the absence 
of external fields | ]28|j . The result is that the value of the induced mass m, obtained by 
keeping the momentum dependence of the gap function inside the pertinent integrals in the 
Schwinger-Dyson equation, is half the value of the mass gap obtained under the "constant 
mass approximation", i.e. the zero-field limit (B — > 0) in ( |5.13|) should be m = g 2 /47r. 
This should be taken into account in quantitative analyses of the phenomenon, and possible 
detailed applications in condensed-matter physics, which, however, go beyond the scope of 
the present work. 

VI. BEYOND THE QUENCHED APPROXIMATION 

To take into account the contribution of internal fermion loops we begin with a study 
of the one-loop vacuum polarization graph in QED3 in the case of even number of fermion 
flavors. The polarization tensor in the one-loop approximation is given by : 

iWp) = -9 2 H0f tr {-rAQjuSik - P )) (6.1) 

where S(k) is the fermion propagator in the presence of the external magnetic field (|4.3|) . 

The calculation of the polarization tensor is straightforward. Due to the fermionic loop, 
the effects of the translational invariance breaking phase factors in (|4.2|) cancel, and one can 
go directly to momentum space as in the zero external field situation. For our case, the 
photon polarization can be obtained easily by performing a dimensional reduction of the 
four- dimensional result [^,^T|]. We end up with 

IV(p) = {p 2 5^ -p^p u )N {p) + {p\5^ - p± fl p± u )N 1 {p) = p 2 P^N {p) +p 2 P±^N 1 {p), 

(6.2) 

where p 2 = p\ with p± = (0,pi,p2) and 
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00 ds f +1 



87T2 J •'-I 

g 2 z" 00 rfs 



sinh z 



S7T2 JO 



<it>e 



-S</>0 



2z 



s J-l 



sinh z 



[cosh zv — v coth z sinh zt> ] 
[cosh z — cosh zu] — N (p) 



(6.3) 



with 



, 2 1 — v 2 2 cos h -2 — cos h ^ 2 
6 o= m + — z — Po + — 77— : — — P ' 



2z sinh z 



(6.4) 



An outline of the derivation of the above-mentioned formulae is provided in the appendix. 

For weak magnetic fields, we will have \J~eB « S(0), where S(0) is the dynamically 
generated fermion mass. Note that it is the opposite to the limit encountered in the case for 
the strong magnetic field |],^(]]. In the weak-field limit, we can expand the above functions in 
a power series of z — s£ 2 (0) ( ^I^) ) an d take the leading and next to leading order behavior 
as z — > 0. 

We have, as z — > 0, the following expansions to order e 2 B 2 : 



1 - v 2 z 2 

Vj2 _|_ 1 u J2 * 1 1 ,:> \2,^2 



4 



p 2_ ^ (1-V 2 ) 2 P 2 + 0(Z 4 



No(p) 



9 f ds f dve- s ^(l-v 2 ) 



Itt2 Jo 



Ni(p) 



fj 



s J-l 

2 rco ds C+ 1 
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1 - £<!-«■) 



+ 0(z 4 



7T2 JO V s J ~ l 



dve 



-s4>o ■ 
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-3 + 2v 2 + v 4 ) + 0(z 4 



(6.5) 
(6.6) 

(6.7) 



Note that when z — > the term Ni(p) vanishes and we recover the usual form for the 
polarization tensor. 

To simplify the integrals, we also expand the exponential e~ s ^° in a power series in z so 
that, 



V 48 



1 - v 2 ) V 



(6.8) 



This simplifies the s-integrals. We then end up with: 



No(p) 



9' 



-[(: 



2E Z 



sin 1 k + 



2itp lx 2 p 2 ' p 

2e 2 B 2 r 2Sk 2 (3 + 2k 2 ), 
■{sin K } 



1 1 3p 

. , 2S/t 2 (15 + 10/t 2 + 8/t 4 )„ 
l sm K — )] 



p 

5e 2 B 2 p 2 



2p { 



15p 



(6.9) 



g 2 e 2 B 2 



Ni(p) 



2i\p 



— [sin k — p- 



3 — k 2 + Qk 4 
~6£ 



(6.10) 
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where k 2 = 4S f + 2 . Note that for physical processes < k < 1. When, eB = and £ = 
(i.e. k = 1 )we recover the known results (equation |3.7|) , namely N (p) = — A^(p) = 0. 



On the other hand, in the presence of the magnetic field one can readily see that in 
the limit £ — > the function Ni(p) blows up when k — > 1 (i.e. the massless case ) due 
to the presence of the factor E in the denominator. However, for a super renormalizable 
theory this seems unphysical. A resolution to this puzzle can be provided by the generation 
of a dynamical fermion mass £ ( however small ) in the presence of the magnetic field. 
This observation points to the magnetic catalysis even in the case for the weak fields in 
three-dimensional QED. Generation of such a mass would prevent the appearance of the 
divergences. 

For E 7^ 0, No(p) and Ni(p) behaves, when E << p and p — ► oo, as: 

Ar/ , g 2 ( Ae 2 B 2 5e 2 B 2 sm 2 9\ . . 



2vrp 5 V 2 3E 

where we have used the parametization k 2 = k 2 sin 2 9, as before. Hereafter we will be using 
the approximations (|6.11| ) and (|6.12 ) which are justified in the weak-field case. 



As in the case with the quenched treatment let us proceed to get the gap equation. Again 
we shall adopt the approximate qualitative approach of the previous section as regards the 
effects of the translational invariance breaking phase factors in (|4.2D due to the presence of 
the external field. To take account of the photon propagator, we can invoke the large-N 
argument to sum up the photon propagator in the ladder approximation : 

DAv) = DIM + D^( P )^p(p)D° p Ap) + ■■■ (6.13) 

where D° is the free photon propagator. To facilitate our calculations let us use the Landau 
gauge for the zeroth order propagator, so that 

DIM = (6.14) 
Using the algebraic properties of the projectors 

p^MPu P ip) = p, p (p) 

P^(p)P±up(P±) = P±up(P±) 

P±,M)P±U P (P±) = P±u P (p±) (6.15) 
we can sum the series in (|6.13|) to get 



D,M ' 



p 2 (l-NM) 



Pfii/ ~l" P± P u 



(1 - (nm + %nm) 



(6.16) 



To go beyond the case of quenched approximation which we discussed in the previous section, 
we need to include the polarization effects in our analysis treatment. To perform this 
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improvement we replace the photon propagator in (|5.3|) by the ladder resummed one, given 
by ( |6.16|) . For the fermion propagator we proceed as we did previously for the quenched 
case: starting from equation ( |5.4|) we expand the term e" 



sk 2i - 



e -sk 2 . e - S fc 2 (iS£h£- 



and get the expression 

Then we expand the second exponential in powers of z = eBs, retain 
the terms which are at most of second order in eB and integrate over s. Taking traces over 
the Dirac 7 matrices we finally get: 



d 3 k 



2e 2 B 2 k 2 



(2n) 3 (k 2 + k 2 + J: 2 (k)) 



1 + 



1 



(E 2 + k 2 ) 3 (p- k) 2 (l - N (p - k)) 



x 



2 + 



(p-fc) 2 



(6.17) 



Let us now set p = and as in the case of the quenched aprroximation let us make the 
substitution £(&;) — S(0) = m to get 



d 3 k 



2e 2 B 2 k 2 
(m 2 +k 2 ) 3 



(2tt) 3 {kl + k 2 + m 2 ) k 2 (l - N (k)) 



x 



2 + 



(NoW + ZN^k)) 



(6.18) 



Now as k 2 = /cq + k 2 and we can use the parametrization k = kcosO and \k\ = ksin9 . Let 
us write, using ( |6.11| ), 

N (p) = E(p) + F(p) sin 2 9 (6.19) 
where E(p) = —^-(1 — 4g2 f 2 ) and F(p) = Accordingly, we can rewrite ( |6.18|) 



as 



dk 



dO- 



sin 9 



(2vr) 2 (k 2 + m 2 ) Jo (1 - E(k) - F(k) sin 



x 



2 + 



iVi(fc) sin 2 fl 



1 - E (k) - (F(k) + iVi(fc)) sin 2 9) 



2e 2 B 2 k 2 sin 7 9 
(m 2 + fc 2 ) 3 



(6.20) 



The angular integral can be performed by making a change of variables y = cos 9, so that 
we end up with: 
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(6.21) 



where 6 2 



and a 2 



l-E-F-N! 
Ni+F ' 
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However, this equation is difficult to handle, so we prefer instead to expand equation 
Q6.20Q in powers of e 2 B 2 before doing the angular integration. After some rather tedious 
algebraic manipulations we end up with: 



eB _ f 2 1-A 

V = ^ = inTYi (6 ' 22) 

The quantities appearing in the right hand side are given by the following expressions: 

2 f°° x 2 
nfi Jo (x 2 + l)[x 2 + h(n, x)] ' 

1287T 3 f°° x 4 

Tl = ^rJo dx (x 2 + irh^xy (6 ' 24) 



T 128tt 3 12x 4 + x 6 

2 V Jo (4 + x 2 ) 3 (x 2 + l)[h(/i,x)} 2 ' 

In the above expressions we have used the notations / = 27r^P and /i = 27Tp, already used 
in the previous section; moreover, we have employed the expression: 



h(ii,x) = - 



x 2 — 4 i x 
1 H sin ' 



2x ' 



Numerical computation of the integrals yields ^ as a function of which can be used 
to produce the lower branch in figure [| This represents the solution of the SD equations in 
the region of small magnetic fields. (Note that only the small- .B part of the curve is depicted). 
We see that the dynamically generated mass in this case is substantially smaller than in the 
previous section. Presumably this reflects the fact that, due to the Pauli principle, the 
condensate tends to decrease. We note that in the quenched case the back reaction of the 
fermions is not really felt, so this fact has no consequences in that case. 



VII. LATTICE FORMULATION 

We now proceed with a description of the lattice formulation of the problem. The lattice 
action is given by the formulae given below. 

S = V E F^(x)F^(x) + *nQn,n^n> (7.1) 

n ,ri 

F )lu (x) = a^(x) + al(x + //) — a^(x + v) — a s v (x) 

Qn,n> = $n,n> ~ K J2i S n> ,n+fi(r + T/O^Vn/l + <*n',n-/i(r - l^U^^V^^}. 
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The indices n, n' consist actually of three integers each, (m, n 2 , n 3 ), labeling the lattice 
sites, while ji denotes directions, r is the Wilson parameter, K the hopping parameter, U n ^ = 
e l 9 aa n(i^ y n ^ = e ^aA njl ^ a s re p resen t s the statistical gauge potential and the external 
electromagnetic potential. (5q = is related to the statistical gauge coupling constant 
in the usual way. On the other hand, we denote by e the dimensionless electromagnetic 
coupling constant of the external electromagnetic field. In our treatment we will use naive 
fermions, so we set r = 0. Initially we will consider a homogeneous magnetic field; thus one 
should construct a lattice version of the homogeneous magnetic field. This has already been 
done before in in connection with the abelian Higgs model. We more or less follow them, 
but follow a slightly different prescription, which we describe below ||. 

Since we would like to impose an external homogeneous magnetic field in the (missing) 
x 3 direction, we choose the external gauge potential in such a way that the plaquettes in 
the x±x 2 plane equal B, while all other plaquettes equal zero. One way in which this can be 
achieved is through the choice: Az(ni, n 2 , n 3 ) = 0,for all ni, n 2 , n 3 , and 

l),m ^ N,A 1 (N,n 2 ,n 3 ) = -f (N + l)(n 2 - 1), (7.2) 
1), n 2 £ N, A 2 (n h N, n 3 ) = +f (N + l)(m - 1). (7.3) 

where N 3 is the number of points on the (cubic) lattice. It is trivial to check out that all 
plaquettes starting at (n 1; n 2 , ra 3 ), with the exception of the one starting at (N, N, n 3 ), equal 
B. The latter plaquette equals (1 - N 2 )B = B - (N 2 B). One may say that the flux is 
homogeneous over the entire x\x 2 cross section of the lattice and equals B. The additional 
flux of —(N 2 B) can be understood by the fact that the lattice is a torus, that is a closed 
surface, and the Maxwell equation V • B = implies that the magnetic flux through the 
lattice should vanish. This means that, if periodic boundary conditions are used for the 
gauge field, the total flux of any configuration should be zero, so the (positive, say) flux B, 
penetrating the majority of the plaquettes, will be accompanied by a compensating negative 
flux —(N 2 B) in a single plaquette. This compensating flux should be "invisible", that is it 
should have no observable physical effects. This is the case if the flux is an integer multiple 
of 27r : N 2 B = m27t — » B = mj^, where m is an integer. Thus we may say (disregarding the 
"invisible" flux) that the magnetic field is homogeneous over the entire cross section of the 
lattice.0 The inte ger m may be chosen to lie in the interval [0, ^-], with the understanding 
that the model with integers m between ^ and iV 2 is equivalent to the model with integers 
taking on the values iV 2 — m, which are among the ones that have already been considered. It 
follows that the magnetic field strength B in lattice units lies between and tt. The physical 
magnetic field B p h ys is related to B through B = ea 2 B p h ys , and the physical field may go to 
infinity letting the lattice spacing a go to zero, while B is kept constant. 

An important remark is that the magnetic field is not allowed to be too big in lattice 
units, since then the perturbative expansion of the expressions e ieaAnf * would yield significant 



3 To check this translational invariance we measured the fermion condensate at every point in the 
x\x 2 plane. The results were the same at all points within the error bars, confirming homogeneity. 
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B 2 , B 3 , . . . contributions with the accompanying vertices, in addition to the desirable terms 
which are linear in B. A trivial estimate of the critical field strength is obtained from the 
demand that the cyclotron radius corresponding to a given magnetic field should not be 
less than (say) two lattice spacings. This trivial calculation yields B < |. Of course the 
above limitations apply strictly only to the case where the statistical gauge field has been 
turned off; in the "interacting" case, one does not really know whether there exists a critical 
magnetic field, after which discretization effects are important. With this remark in mind, 
we depict in the figures of the following sections the results for the whole range of the 
magnetic field, from to ir. 

For the fermion fields we used antiperiodic boundary conditions in the time direction 
and "fixed" boundary conditions in the spatial directions; the latter boundary conditions 
mean that we consider fermion fields vanishing on the boundaries. 



VIII. LATTICE RESULTS 

A. Zero Temperature Results 

We are now going to present the results pertaining to the T = case. The first set consists 
of measurements of the fermion condensate versus the magnetic field for a 16 3 lattice in the 
strong coupling regime for the statistical gauge field (/3g = 0.10) for three values of the bare 
mass (figure |j). Before going on with the specific features of this results, let us remark that 
to facilitate comparison with the analytic results we measured the magnetic field in units 
of its maximal value: thus we used the parameter b, defined by: b = R B = e - B p fe ^° — . 

i3max €.^phys^ \ max 

Since B max = n, as explained previously, we get: b = — and b runs from to 1. We see in 
figure f| that for all three masses the plot consists of two parts with qualitatively different 
behaviour. For b smaller than about 0.3 we find a dependence of the condensate on the 
external magnetic field, which is nearly linear, however, in view of the analytical results 
obtained in section 5 about the quenched case, we understand that we see the quadratic 
behaviour found there; however, there is a negative quartic term coming into the game, as 
we also saw, and this "straightens out" the quadratic curve and makes it almost linear. For 
big magnetic fields we find points that could possibly be fitted to a logarithmic type of curve. 
The logarithmic dependence 



a 



a 



9 2 



has been found JT^| by an approximate solution of the Schwinger-Dyson equations in the 
regime of strong magnetic fields. We have included such a logarithmic fit for m = 0.050 
in figure f|. In addition, for this mass some points in the intermediate region are included. 
They show a smooth interpolation between the two regions. Thus in both the strong and 
weak magnetic field regimes we find a nice qualitative agreement of the analytical solutions 
with the Monte Carlo results. In the figure we have also included the extrapolation to the 
zero mass limit. 

The magnetic field has been characterized as "strong" or "weak" through its comparison 

2 

with the fine structure constant a = f- and the dynamically generated mass m ~ S(0). 
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Since figure £| contains the strong coupling data, it would be interesting to explore the fate 
of the whole picture shown in figure ^ as the gauge coupling g moves away from the strong 
coupling regime. One would naively expect that the magnetic fields will be more easily 
characterized as "strong" , as compared to the smaller and smaller scale set by the gradually 
weaker coupling constant. Thus, the almost linear part should be restricted to the very 
small magnetic fields and eventually disappear. This is what one may see in figure [5], which 
is similar to figure [|, the only difference being that the gauge coupling constant is now in 
the intermediate coupling regime, rather than the strong coupling of figure |j. We see that 
the almost linear part is now restricted in the region between b = and b ~ 0.12 — 0.15. We 
may also have a semi- quantitative estimate of the new "critical" magnetic field b c , defined 
as the maximum b which fits into the almost linear behaviour. Inspired by the inequality 
eB « (f^) 2 , let us suppose that eB c = x(f^) 2 , with x a very small number; we also make 
the further assumption that x does not depend on g. We will check crudely whether this 
assumption is reasonable given our results. Converting everything to lattice units, we find 
that b c = j^sjr- From this we infer b C2 = b Cl (^) 2 . Using (3 Gl = 0.10, [3 G2 = 0-15 and 
b Cl = 0.3, we find for b C2 , the value 0.13, which is surprisingly close to the value given by the 
data of figure E. Of course, the 4r dependence of b c tells us that the weak field behaviour 

P G 

will be even more suppressed as we move towards the weak gauge coupling; this is what we 
have seen in direct simulations in this regime. Thus, it is plausible that the dependence of 
the "critical" magnetic field has a 4r dependence. 

We now make contact with the results of ||, where we studied the model with the 
statistical gauge field turned off. We had found there that for big enough b the condensate 
stopped showing a monotonous increase with b, at b = 0.5 it had a local minimum and then 
had a succession of maxima and minima, up to b = 1. Moreover, there was a spectacular 
volume dependence. One expects, of course that this "free" case will be reached for big 
enough (3 G . In figure |6| we show the results for /3 G = 0.5 and (3 G = 1-0 for various volumes. 
For (3 G = 0.5 the curve shows the first sign of non monotonous behaviour at b = 0.5, while 
at f3 G = 1.0 the succession of maxima and minima is clear. However, there is no detectable 
volume dependence, so we can be sure that, even at this large (3 G) the limit of switching the 
gauge field off has not yet been reached; it will presumably be reached for even bigger values 
of (3 G . One should add that in the "free" case the role of the bare mass is very important, 
since it is eventually the only source of mass generation. This is at the root of the large 
volume dependence showing up in the "free" case: at fixed volume the condensate goes over 
to zero for vanishing bare mass. In the full model, though, the interaction with the gauge 
field generates a dynamical mass, independently from the value of the bare mass. This is why 
in the "interacting" case the volume dependence is small, permitting a smooth transition to 
the thermodynamic, as well as to the massless, limit. 

The simulations are done at finite values for the (bare) mass; the massless limit is taken 
by extrapolating the results for several bare masses to the limit m — > 0. Figure [7| shows the 
process of this extrapolation for three values of the gauge coupling constant. The external 
magnetic field has been set equal to a typical value {b = 0.188); the picture is similar for 
all values of the magnetic field strength. For f3 G = 0.10, which lies in the strong coupling 
region, the extrapolation is linear with negative slope. This line is pointing to a relatively 
big value for the condensate in the chiral limit. For somewhat weaker coupling (/3 G = 0.20), 
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the curve is still a straight line, but the slope is positive and it points to a smaller value at 
m — > 0. In both of these cases the mass dependence is not very pronounced, because it is 
the strong gauge coupling which dominates in the formation of the condensate. In the third 
case (Pg = 0.35), which lies in the weak coupling, one can no longer fit to a straight line; a 
quadratic fit has proved necessary for all /3q values smaller than 0.2. 

Figure |8] contains the zero mass limit of the condensate (obtained through the procedure 
illustrated in figure [7]) versus Pg, for four values of the external field. We observe that in the 
strong coupling region the b-dependence of the condensate is rather weak; on the contrary, 
at weak coupling, the external field is the main generator of the condensate, and we find an 
increasingly big b-dependence, as we move to large Pg- Note that the biggest value of b we 
have used in systematic measurements, such as the ones in figure |§, is 0.3. This takes into 
account that for larger values of b the function of the condensate stops being monotonous for 
large /3g, as may be seen on figure |. Thus we have restricted our study to a b region which 
is safe for all values of the coupling. From this preliminary quenched study we conclude that 
a non- vanishing value for < > develops for this small volume even at weak coupling in 
the presence of an external magnetic field. We have not tried to study systematically the 
approach to the continuum limit. 

B. Finite Temperature Results 

We expect that the fermion condensate, generated by any mechanism (explicit mass 
term, gauge interactions, external fields) should vanish at high temperatures. This is the 
study we now turn to: we employ asymmetric lattices and consider the behaviour of the 
condensate versus Pg- Before describing the behaviour of the condensate, let us first see the 
Pg dependence of the Wilson line. 

Figure |9] depicts the Wilson line versus Pg for lattices of temporal size N T = 4 and 
various spatial volumes. We observe in the figure not only the decrease of the Wilson 
line with increasing spatial volume, but also the fact that initially this quantity is almost 
independent of Pg, but at Pg — 0.25 its dependence on p G starts showing up. It is important 
that this value of Pg is independent from the spatial volume. 

In figure we show the relationship between Wilson lines on lattices with N T = 2 versus 
N T — 4 and N T = 6. The result for the 16 2 * 6 lattice lies "below" the result for the 16 2 * 4. 
This is due to the fact that the former lattice is closer to the zero temperature (symmetric) 
lattices, as compared to the latter. The value of Pg, above which there is Pg dependence is 
substantially bigger for iV T = 6 than for N T = 4 or iV T = 2. The Wilson line for the iV T = 2 
case approaches an asymptotic value for large Pg- This is not very obvious in the other two 
cases, because they lie farther from the infinite temperature limit. Also in this case the 
statistical fluctuations are very large, resulting in big errors. With this in mind we have put 
in the figure only the errors for the case N T = 2.. 

Figure [IT] contains the zero (bare) mass extrapolations for the condensate as a function 
of Pg- The external magnetic field is set to b = 0.1. The uppermost curve contains the 
results for a symmetric lattice (16 3 .) It is easily seen that it is a quite smooth curve and 
presents no apparent discontinuities of any sort. The data for the asymmetric lattice 16 2 *4 
follow the ones of the symmetric lattice at strong coupling; in the weak coupling region the 
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condensate for the asymmetric lattice appears substantially smaller than its 16 3 counterpart. 
This is what one should explain on account of the symmetry restoration scenaria at finite 
temperatures. The 16 2 *4 data can be described by two branches, one containing the strong 
and the other the weak coupling results; the two branches join at about j3c = 0.4, but their 
slopes do not coincide. There is a discontinuity at this value of /3g, which we interpret as 
the symmetry restoring transition at finite temperature. On the same figure we have put 
the results for a lattice of bigger spatial volume (24 2 * 4), away from the "critical" (3g = 0.4 
value. These data do not differ substantially from the ones for 16 2 * 4. 

Figure ^ contains results similar to the ones of figure |TT], but the value of the external 
magnetic field differs: b = 0.305. The same basic picture appears here, as well: we may again 
spot the discontinuity at (5q = 0.4 for N t = 4. In addition to the data of figure [TT] we have 
put the data for a 24 2 * 6 lattice, which is expected to lie closer to zero temperature. The 
data are smoother than the ones for 16 2 * 4 and they lie much closer to the 16 3 results; this 
makes it more difficult than before to spot a sudden change in slope; however, this change 
is present even in this case. The new element here is the data for N t = 2. The change in 
slope here is very pronounced and substantiates our claim that we have a phase transition 
around (3q = 0.4. 

Since we now have data for several asymmetric lattices, we are in a position to show the 
temperature dependence of the condensate. This is done in figure 13, in the weak coupling 
regime, for two values of the magnetic field. The zero mass extrapolation of the results has 
been used and the temperature in lattice units is as usual. We observe the fall of the 
condensate at high temperatures, which is more dramatic for the smallest value of b. This 
figure is of the same qualitative form with figure 2 of || , which was derived analytically for 
the case where we had no statistical gauge field at all. 

In figure |T4| we show the time evolution of the condensate for two values of the magnetic 
field at weak coupling for a symmetric lattice. The important feature here is the very small 
magnitude of the statistical fluctuations, resulting in relatively small errors. 

The situation in figure |14] changes in the asymmetric lattices and the results of figure |14] 
should be contrasted against the ones in figures 15 and 16. In figure 15 we have exactly the 
same parameters as in figure ^4[ however now we have a 16 2 * 4 lattice. It is evident that 
the fluctuations have grown about one order of magnitude larger. 

The statistical fluctuations grow even larger for the 16 2 *2 lattice, whose results are shown 
in figure [16| Moreover, this figure gives a feeling of the way the average of the condensate 
is approaching zero at high enough temperatures. The outcome spends most of its time at 
small values and has some exceptional big spikes from time to time; these latter become 
more and more rare as the spatial volume increases. 



C. Non-uniform Magnetic Field 

In the previous sections we have considered the case of a uniform external magnetic 
field. There is however potential physical interest in the effects of non-uniform fields, which 
become important in case the above model has relevance to the physics of high-temperature 
superconductors. Indeed, it will be of interest to examine the effect of electromagnetic vortex 
lines on the induced fermion (holon) condensate at the nodes of (i-wave superconducting 
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gaps . A simple-minded model for such vortex lines could be that of flux tubes of magnetic 
field. The full problem would be to take into account interactions among the flux tubes, 
which could shed light also in the confining aspects of the gauge groups in three dimensions. 
It is only by lattice methods that one may treat the problem, in view of the very big 
computational difficulties in the analytical approach. In this first treatment of the problem 
we switch off the fluctuating statistical gauge field and consider the response of the fermions 
to the background field only. A full treatment of the problem, including the statistical gauge 
interactions is left for the future. 

Let us describe the technical procedure to construct a non-uniform magnetic field on the 
lattice. We will consider M x M plaquettes parallel to the x±x 2 plane, around the center 
of the lattice, which will be penetrated by magnetic flux equal to B each. The remaining 
plaquettes will not carry any flux. Then we are going to measure the condensate at the 
center (y, y, y) and along a line passing through it and consisting of the sites (—, y + 1, y), 

(N N_ 1 r> N_\ (N at N_\ 

\ 2 ' 2 2 )■>■■■ ■> \ 2 ' ' 2 /" 

The fact remains that the total flux through the lattice should be zero, because of 
divli = 0. Thus, for each flux B penetrating a given plaquette, there should be an opposite 
flux somewhere else in the lattice. To construct the magnetic field that we mentioned above, 
we followed the strategy to build it up plaquette by plaquette taking care that we put 
the compensating opposite flux through the plaquette starting at the point (N,N,n 3 ). If 
we already have a configuration of gauge potentials on a lattice, the procedure to add a 
plaquette of flux B at the plaquette at (iVi, N 2 , n 3 ) (with the corresponding compensating 
flux — B at (N, N, n 3 ) ), consists of adding to the preexisting links the quantities denoted by 
AA k below. AA 3 (rii,n 2 ,n 3 ) is set to zero for all values of the integers ni,n 2 ,n 3 . AAi and 
AA 2 are also set to zero, except for the links where an explicit different statement is made. 
For the plaquette starting at the site (N 1: N 2 , n 3 ), Ni 7^ N, N 2 7^ N, we choose: 

AA 1 {N 1 ,n 2 ,n 3 ) = -B, n 2 = N 2 + 1, N 2 + 2, . . . , JV, n 3 = l,...,N 

AA 2 ( ni ,N 2 ,n 3 ) = +B, n 1 = N 1 + l,N 1 + 2,...,N, n 3 = 1, . . . , N 
For N\ = N we impose: 

AA 1 (N,n 2l n 3 ) = -B, n 2 = N 2 + 1, N 2 + 2, . . . , N, n 3 = l,...,N 

AA 2 (ni,n 2 ,n 3 ) = everywhere, 

while for N 2 = N : 

AAi(rii, N, n 3 ) = everywhere, 

AA 2 (ni,N,n 3 ) = B, n x = i\q + 1, Ni + 2, . . . , JV, n 3 = l,...,iV 

In the following we consider the model with the statistical gauge field turned off. We start 
with vanishing gauge potentials everywhere on the lattice, go through the plaquettes in the 
central region and add the above AA^ quantities to the corresponding links. In this way 
we end up with the flux B in the central plaquettes and the compensating flux for all the 
plaquettes at (N,N,n 3 ). The flux through these latter plaquettes should be "invisible", as 
explained in section , so B must take the values jrsn, n = 0, 1, . . . , 
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In figure [17] we show the results for a central region of non- vanishing flux of extent 6x6. 
More specifically, for the 16 3 lattice we have been using, the region with constant non-zero 
flux contains the plaquettes starting at (rii, 712, ^3), with 6 < n\ < 11 and 6 < n 2 < 11, 
while ns takes all values. Note that nothing depends on the value of n 3 . The uppermost 
curve in the figure depicts the result for the condensate at the site (9,9,9). The remaining 
curves represent the corresponding results for the sites (9, 12,9), and (9, 16,9). The curves 
with ri2 = 10, ri2 = 11 are quite similar to the 77-2 = 9 curve. The first substantial change 
takes place at the site (9,12,9), which lies exactly on the boundary of the above region. 
The remaining curves dive together to a value which is accounted for by the explicit mass 
term and has very little to do with the external magnetic field. Thus, we find a drop in the 
condensate value taking place exactly on the boundary of the central region. 

To visualize the fall of the condensate on the boundary, we fixed the b parameter to 0.111 
(a typical value) and plotted the value of the condensate along a straight line passing from 
the center of the lattice. We find the symmetric bell-shaped plot shown in figure 



IX. CONCLUSIONS 

In this work we have studied in detail, by means of analytic and lattice methods, the 
phenomenon of magnetic catalysis in even-flavour QED 3 , namely the magnetic induction of 
a chiral-symmetry breaking fermion condensate as a result of the influence of an external 
magnetic field. We have shown that the scaling behaviour of the induced condensate with the 
external field varies according to the strength of the latter. In the weak-field regime, there is 
a quadratic increase of the condensate with increasing external field, to be contrasted with 
the logarithmic scaling behaviour in the regime of strong external magnetic fields. However, 
it seems that the transition from weak to strong fields is smooth, at least as far as the 
induced condensate is concerned, and we would characterise it as a cross-over rather than 
a phase transition at some critical value of the external field. This constitutes a prediction 
of the gauge theory, and it may be tested in experiments of relevance to high-temperature 
superconducting materials. It would be interesting to repeat the (lattice) computations for 
the case of four-fermion contact interactions to check on this behaviour. This would differ- 
entiate between the two models as possible candidates for the nodal spin-charge excitations 
in d-wave high-temperature superconductors |7|.|T5l . 



It should be stressed that the analytic methods that lead to this scaling are approximate, 
and should be considered only as giving a qualitative treatment of the phenomenon. The 
important effect of the external field is the breaking of translational invariance on the spatial 
plane, and this leads to technical complications in solving the pertinent Schwinger- Dyson 
equations in the case of weak magnetic fields, mixing configuration and momentum space 
integrals. To bypass this problem, we adopted an ansatz for the fermion propagator in 
the presence of a weak external field, which although maintains formally a translational 
invariant look (in the sense of its being expressed in terms of a Fourier transform depending 
on a single momentum variable), however it incorporates the effects of the magnetic field 
in the pertinent coefficients. Comparison with the quenched lattice results showed that 
the predicted scaling of the induced condensate with the magnetic field is (qualitatively) 
captured by this ansatz. 
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In addition to the uniform external field case, we have also presented preliminary 
quenched lattice results in the case of flux tubes of magnetic field. This situation might 
also be of relevance to realistic situations in high-temperature superconductors, as being 
related to the effects of electromagnetic vortex lines on the opening of a fermion gap at 
the nodes of the superconductor, within the context of the gauge theory approach [7|]. Our 
results in the non-uniform magnetic field case have indicated that the fermion chiral con- 
densate is non zero and scales with the magnetic field of the flux tube at the core of the 
latter, but decays very fast outside the tube. Our considerations did not properly take 
into account interactions among flux tubes. The latter is an important issue, which might 
also bear some relation with the issue of confinement of the three-dimensional theory. We 
expect that a proper treatment of this problem will become available only upon the use of 
dynamical fermions on the lattice. 

Another important issue we would like to address for future work is the computation 
of thermal conductivities in the context of the model of section 2, used in our simulation 



of the physics of planar high-temperature superconductors. As discussed in [[T4| there are 
scaling differences of the thermal conductivity between the gauge (QED3) and four- fermion 
models, which would be important to analyse in detail in the context discussed in this work 



for comparison with experiments of high-temperature superconductors JT3|. At present, the 
analysis of the thermal conductivity has been performed in the real time formalism [15,14], 
and the extension to a lattice analysis is not trivial. We hope to return to this important 
issue in a future publication. 
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APPENDIX A: CALCULATION OF THE ONE-LOOP VACUUM 
POLARIZATION IN THE PRESENCE OF THE EXTERNAL MAGNETIC FIELD 

The photon polarization tensor for QED in the presence of the electromagnetic field was 
first performed by Tsai [3(],31]. Here, for the sake of completeness, we outline his calculations 



but we will work in three-dimensions instead of four- dimensions. 

Let us begin with the one- loop vacuum polarization graph (fj). The polarization tensor 
to this approximation is given by 

/d k r ~ ~ 1 

— tr [YS(k)YS(p - k)} (Al) 
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where the Fermion propagator is the one in the presence of the constant external magnetic 

field ra 



S F (k) = i / dse- s(fc o+ m2+fc2 ^H)[(m - 7 • fe) - 2(71*2 - 7^) tanhz](l - i 7l72 tanhz) 
Jo 

= % I" rf S e- s(m2+fc o+ fc2 ^)((m - 7o Jfe )(l - z 7l72 tanh^) - (7 • k)—L-) (A2) 
Jo cosh z 



where z = e.Bs. Accordingly, ( |A2j ) leads us to 

n^(p) = ^ jT d Sae -cxo(«.*)+w.(«*-*)) 



x tr 



1 



7 /x ((m - 7ofco)(l - i 7 i 7 2tanh2;i) - (7 • fe) — -3 — ) 

cosh 

7"((^-7o(p- fe)o)(l - «7i72tanh22) - (7 • (p - fe)) -5—] 

cosh z 2 



where = eBst, i — 1,2 and 



Xo(s, fe) = s(m 2 + feg + fe' 



,tanh z N 



Let us now make the change of variables 

1-v 



si 



1 + v 

S2 = — : — S, 



2 7 2 
with s e [0, 00) and u G [—1, 1]. Accordingly, we get 

Xo(si, fe) + Xo(s 2 ,P - fe) = s[0o(p) + 01 (p, 

with 



(A3) 

(A4) 

(A5) 
(A6) 



0o(p)= ™ 2 + 
0i (p, k)= (fe - 



1 - v 



cosh zt> — cosh z 
2z sinh z 



1 + v n9 tanh zi + tanh z 2 
-^Pof + 



k- 



tanh Z2 



tanh Z\ + tanh z 2 



P 



(A7) 



The loop integrals can then be performed very easily via standard Gaussian integrations 
some of which are listed below: 

d 3 k 



d 3 k 



(2k) 



-s<j>i(p,k) 



(4vr S )2 sinh, 



cosh zi cosh z 2 



e- s ^k = ^p Io 



(2tt) 3 

J (27T)3 e ^ 



(2, 

d 3 k 
(27p 

(2tt) 3 * j 



tanhz 2 
tanh^ + tanh , 2 PJ ° 
1 + f tanh z 2 
2 tanh zi + tanh z 2 



^oPoP 



tanh z\ 
tanh + tanh z 2 



PiPj 



s(tanh z\ + tanh z 2 ) 



(A8) 
(A9) 
(A10) 
(AH) 
(A12) 
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We will also need the following identity for the gamma matrices: 

tr [(1 - r^tanh 2)7^7,,]= -4(5^ - tanhz(<^5£ - 5*8%)) 

= -4(V-tanhz%) (A13) 

where F^ v is the covariant representation of the external magnetic field strength. The other 
traces can also be evaluated by the use of the Dirac algebra ( j2.6|) . 

Putting everything together 0,^l|] one gets ( after an integration by parts ) 

where 

V = [(<W 2 - P»Pv) R o(p) + (8iuP± - P±vP±v)Ri(p)} (A15) 

with 

Ro{p) = (cosh zv — v coth z sinh zv) (A16) 
2 

R\{p) — q— [cosh z — cosli^f] — Ro(p) (A17) 

sinh z 

Note that, unlike its four- dimensional counterpart, the vacuum polarization tensor in 
three dimensions is not divergent and there is no need add any counterterms. One can see 
this by checking the absence of poles at s — > 0, which is the place where poles usually show 
up in proper time methods. 
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FIGURES 




FIG. 1. One- loop vacuum polarization for photons (wavy lines) in QED%. The solid lines with 
crosses represent fermions in the presence of an external magnetic field. 




FIG. 2. The Schwinger-Dyson equation for the fermion self-energy. The curly line indicates 
the Ug(l) statistical photon. Solid lines with crosses represent fermions in the presence of the 
external magnetic field. Blobs indicate quantum corrections (loops), which are ignored in the 
ladder approximation. Quantum dynamics of the electromagnetic field has been suppressed. 
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Mass Gap from numerical solutions of SD equation 
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FIG. 3. Solution of Schwinger-Dyson equations for the quenched and dynamical fermions. 
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fermion condensate for m=0.1 
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fermion condensate for b=0.188, vol=16 A 3 
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Wilson Line for 3d non-compact QED 
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fermion condensate versus T (zero fermion mass, betag=0.60, L_s=16) 
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FIG. 13. < tyty > versus the temperature for two values of the magnetic field. 
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fermion condensate for fermion mass=0.025, betag=0.60, vol=16 A 3 
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FIG. 14. Time evolution of < tyty > for a symmetric (16 3 ) lattice and two values of the magnetic 



field. 
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fermion condensate (bare fermion mass=0.025, betag=0.60, b=0.305) 
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FIG. 15. Time evolution of < Hty > for a 16 2 * 4 lattice. 



fermion condensate (bare fermion mass=0.025, betag=0.60, b=0.305) 




FIG. 16. Time evolution of < > for a 16 2 * 2 lattice. 
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16 A 3, m=0.05, flux only in the central 6x6 area 
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FIG. 17. < > versus magnetic field strength where the flux is non zero only in a central 
region extending over 6x6 plaquettes. The condensate at sites labeled 9, 12, 16 (see text) is shown. 
The corresponding distances from the center of the flux tube are 0, 3, 7. 
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FIG. 18. < > along a straight line passing from the center of the lattice if the magnetic 
field parameter b is set to 0.111. The central region of non-zero flux is 6x6. 
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